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Abstract 

In the framework of a final state interaction model, we show that the so- 
called radial flow, i.e. the almost linear increase of the inverse slope T with 
the mass of the produced particle, is already contained in the initial condition 
- with a slope < > (the so-called strength of the average radial transverse 
flow) which is larger than the measured one. While the precise value of the 
slope depends on the details of the model, the above result has a very general 
basis - namely the increase with increasing pt of the fixed px suppression, in 
the low pt region. 
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Flows of different types are considered as providing strong support to hydro- 
dynamic models and, more generally, as signals of the production of a system in 
thermal equilibrium in heavy ion collisions. In particular the agreement of the hy- 
drodynamic model predictions on elliptic flow with the experimental data at RHIC 
energies constitutes the basis of the claim that a perfect fluid has been produced 
at RHIC [1]. This interpretation has been challenged in [2] where the authors 
claim that the RHIC data on elliptic flow are inconsistent within the hydrodynamic 
model predictions for a perfect fluid (i.e. with no viscosity). Moreover, an alter- 
native description of elliptic flow at RHIC has been proposed [3] in the framework 
of a final state interaction model - similar to the one previously introduced to de- 
scribe strangeness enhancement [4], J/ip suppression [5] and fixed px suppression 
(jet quenching) [6]. The parameters of the model (j>r-shift due to the interaction 
of the produced particle with the medium and the cross-section of this interaction) 
are determined from the data on fixed pr-suppression. Elliptic flow is then due to 
the different absorption of the produced particle for different values of its azimuthal 
angle - due to the different path length inside the overlap region of the colliding 
nuclei, associated to each azimuthal angle. 

The purpose of the present paper is to examine in a critical way the so-called 
radial flow, i.e. the increase of the inverse slope T with the mass of the produced 
particle - which has been observed in nucleus-nucleus collisions both at CERN and 
RHIC [7] energies and is absent in pp collisions. RHIC results for identified particles : 
ir 1 * 1 , K ± , p and p are presented in Fig. 1 in different centrality bins [7]. The increase 
of T with the mass was derived in hydrodynamic models [8] where : 

T = T + m<u t > 2 . (1) 

T is the freeze-out temperature and < u t > 2 measures the strength of the radial 
transverse flow. The fits [7] to the data using eq. (1) are given by the dashed lines 
in Fig. 1. This fit provides support to the hydrodynamics picture. In particular 
the data are consistent with a value of the freeze out temperature To independent of 
centrality. However, although the slope < Ut > 2 increases with centrality, it is worth 
noticing that the most peripheral bin (60 - 92 %) corresponds to an average number 
of participating nucleons of 14. Thus, according to the hydrodynamic interpretation 
a perfect fluid is already produced at this comparatively low density. This naturally 
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leads us to ask whether the origin of eq. (1) is a genuine final state collective effect 
(which should be already present for N part = 14) or it is contained in the initial 
state (Cronin effect, shadowing, etc.). 

The latter is supported by the following remark. Apart from the increase of 
the inverse slope T with the mass described by eq. (1), we see in the data that 
T increases with centrality. This increase is larger when the mass increases. Now, 
the RHIC data [9, 10] on fixed pr suppression indicate that, at fixed centrality, the 
suppression increases with px in the low pt region (j>t < 2 GeV) we are interested 
in this paper. Indeed, it is widely recognized that, due to the combined effects of 
decreasing shadowing and of the Cronin effect, the initial value of the ratio 

Rau au = dN Au Au /dy(b, p T )/N coU dN™/dy (2) 

in absence of final state interactions increases rapidly with pt- At pr ~ 2 GeV all 
calculations of the Cronin effect lead to an initial value of Ra u Au close to unity 
[11]. On the contrary, the final value of Ra u Au at px = 2 GeV is substantially lower 
than unity [9, 10]. Thus, there is a sizable suppression at px = 2 GeV, whereas at 
Pt = the data agree with model calculations [3] without final state interaction - 
indicating that there is no suppression at pr = 0. 

If this suppression results from the interaction with the medium, it produces a 
decrease of T at each centrality. Obviously this decrease will be larger for increasing 
centrality. Therefore, the final state interaction with the medium has an effect which 
is exactly the opposite of the one seen in the data of Fig. 1 - providing a strong 
indication that the observed phenomenon is already present in the initial condition. 

The remaining part of this paper is devoted to a quantitative study of this point. 
The final result is given in Fig. 1 where the full signs and full lines show our results 
for the "radial flow" in the initial condition. 

In refs. [3, 6] the fixed pt suppression of a pion has been described in the 
framework of a final state interaction model. Due to its interaction with the medium 
the pions suffer a pr-shift, 5pr- This produces a loss at a given p T - There is also 
a gain due to pions originally produced at p T + 5p T - Due to the fall-off of the 
Pt distribution, the loss is larger than the gain, resulting in a net pion suppression. 
The loss and gain differential equation can be solved analytically. The solution gives 
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the survival probability of the pion [3]. The same equation describes the survival 
probability of a particle of type i 

Si(y,Pt,b) = exp y) 

Here Ni(pr) denotes the pr-distribution of the inclusive production of particle i. 
With 5pr — > oo, the quantity inside the square bracket is equal to unity and there 
is no gain term. With 5px — the loss and gain terms are identical and Si = 1. 
p(b, y) is the density of the medium : p(b, y) = dN AA / dy(b) / G{b) , where dN AA / dy(b) 
is computed in the Dual Parton Model [12], including shadowing corrections, and 
G(b) is an effective transverse area computed using Wood-Saxon profile functions. 
All details can be found in [3]. p pp (y) denotes the pp density. Within conditions 
specified in [3] the argument of the logarithm is equal to the ratio of freeze-out time 
to initial time. As stated above we shall restrict ourselves to the low p T range used 
in the determination of the inverse slope (p? < 2 GeV). 

The pion suppression data can be described [3, 6] using a = 1.3 mb and SpT = 
p^ 5 /(20 GeV 1 / 2 ). In principle a can depend on the type of particle produced. How- 
ever, as it has been argued in [3, 6], the final state interactions described by eq. (3) 
take place largely at the partonic level - in which case all cr's should have similar 
values. Moreover, it has been shown in the second paper of [3] that the mass depen- 
dence of i>2 for 7r, K and p can be described by assuming that the cross-section for 
the three types of particles are the same. (The decrease of v 2 with the mass is due to 
the fact that the py-dependence gets broader with increasing mass. In this way the 
size of the gain term, relative to the loss one, increases. This produces a decrease of 
both the suppression and v 2 with increasing mass). Thus we shall use in the follow- 
ing the same value a — 1.3 mb for the three types of particles. It is, however, worth 
noticing that if the value of a increases with the mass, the mass dependence of the 
slope < u t > 2 in eq. (1) in the initial condition, resulting from our calculation, will 
increase (see below). Thus, by taking the same a for pions, kaons and protons we 
minimize the strength of the "radial flow" present in the initial condition. 

Let us now turn to the only remaining unknown quantity in eq. (3) namely 
the one inside the square bracket. It is well known that the pt dependence of the 
inclusive spectra in the low px region can be described by an exponential in the 
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variable my — m (where my is the transverse mass) [7] 

d 2 N = 1 ( m T -m 

27rm T dm T dy { } 2nT{b){T(b) + m )) U 6Xp ^ T(b) 

where T(b) is the inverse slope. The results of the PHENIX collaboration for iden- 
tified n ± , K ± , p and p in three rapidity intervals are shown in Fig. 2. The fits using 
eq. (4) are given in Fig. 3 of ref. [7]. The resulting values of T(b) are also given in 
[7] and shown in Fig. 1 of the present paper. 

All considerations in this paper, as well as the experimental results in Figs. 1 and 
2 refer to the mid-rapidity region \y\ < 0.35. We use eq. (4) in order to determine 
the quantity in the square brackets in eq. (3). In this way all quantities in the 
survival probability are determined and we can proceed to the calculation of the 
initial px-distribution before the final state interaction has been switched on. These 
initial distributions are simply obtained by dividing, for each value of pr and b, the 
experimental result by the corresponding survival probability eq. (3), i.e. 

27rmTdmTdy^ ^ 2irm T dm T dy^ ^ * 2/'^ T ' ) ( 

The results for the initial pt distribution obtained in this way are shown in 
Fig. 2 (full lines). The resulting values of the initial inverse slope T Initml for the 
various types of particle, in three centrality bins, is given in Fig. 1. We see that 
the inverse slopes of the initial p^-distributions are larger than the measured ones 
and their differences increase with centrality. We also see that the "radial flow", i.e. 
the increase of T with mass is already present in the initial condition. Actually, the 
slope < u t > 2 of the mass dependence, eq. (1), is larger in the initial condition than 
in the data. 

Note that in the calculations presented above we have computed the quantity 
inside brackets in eq. (3) using expression (4) with the measured values of T. It 
is easy to check that using instead the values T Imtial , corresponding to the initial 
Pt distributions, the results are quite similar. Actually, one should solve the gain 
and loss differential equations numerically (instead of using their analytic solution, 
eq. (3)). In this way the values of T would evolve from T Initial at initial time To 
to the measured values of T at freeze-out time 77. Obviously, the result will be 
intermediate between the ones obtained using the analytic solution, eq. (3), with 
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either the initial or the final values of T. Since the difference between these two 
results is only of a few percent it is not worth solving the equation numerically. 

Our consideration have been restricted to light particles n"* 1 , K ± , p and p. How- 
ever, it is natural to expect that the increase of T imtial with mass will saturate for 
larger masses. Indeed, the shadowing corrections decrease with increasing mass and 
the Cronin effect is expected to saturate for large masses. In our approach if such a 
saturation exists in the initial condition it will also exist in the final one, i.e. after 
final state interaction. Thus, particles such as f2, J/ip, ... are expected to show no 
"radial flow" while we expect for them a non- vanishing value of [3]. 

Another obvious prediction of our approach is a substantial "radial flow" in dAu 
collisions. 

In conclusion, experimental data on Au Au collisions at \fs = 200 GeV (as well 
as SPS energy ones) show that the inverse slope T increases with centrality and this 
increase is larger with increasing mass of the trigger. In a final state interaction 
model that reproduces the observed fixed px suppression, the final state interaction 
produces a decrease of the inverse slope and this decrease is larger with increasing 
mass of the trigger. This indicates that the so-called "radial flow" is an effect which 
is present in the initial condition - rather than a collective phenomenon resulting 
from the interaction with the hot medium. 
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Figure 1: Inverse slope T versus mass for ir + , K + and p (left) and ir~ , K~ and p 
(right) in three centrality bins for Au Au collisions at y/s = 200 GeV/c. The data 
points (empty signs) and the dashed lines are from ref. [7]. The full signs and full 
lines are the corresponding values for the initial condition, calculated from eq. (5). 
The circles correspond to central (0-5 %), the squares to mid-central (40-50 %) and 
the triangles to peripheral (60-92 %) collisions. 
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Figure 2: Final (data points) and initial (lines) transverse mass distributions for 
7r + , K + and p (left) and 7r~, K~ and p (right) in three centrality bins for Au Au 
collisions at y/s = 200 GeV. The data are from [7]. The lines are calculated using 
eq. (5) and correspond to the initial condition. As expected, the effect of the final 
state interaction is important for central collisions and negligibly small for the most 
peripheral ones. 
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